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In this work we make a systematic analysis of the correlated processes J/ip — > 77/(1440)//i (1420) 
with 7/(1440)//i(1420) — > KKn, nnn and 3tt, where the role played by the so-called "triangle 
singularity mechanism" (TSM) is clarified. Our results agree well with the experimental data and 
suggest a small fraction of /i(1420) contributions in these processes. This study confirms our 
conclusion in [Phys. Rev. Lett. 108, 081803 (2012)] that the dynamic feature of the TSM can be 
recognized by the strong narrow peak observed in the nn invariant mass spectrum of 7/(1440) — > 37r 
with anomalously large isospin violations. Nevertheless, we explicitly demonstrate that the TSM 
can produce obvious peak position shifts for the same 7/(1440) or /i(1420) state in different decay 
channels. This is a strong evidence that the 7/(1405) and 7/(1475) are actually the same state, i.e. 
£f~) ' 77(1440). We also make an analysis of the radiative decays of 7/(1440) — > (V = <j), p° or uj) 

which shows that such a one-state prescription seems not to have a conflict with the so-far existing 
f*"^ ■ experimental data. Our analysis may shed a light on the long-standing puzzling question on the 

PsJ 1 nature of 7/(1405) and 7/(1475). 
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PACS numbers: 13.75.Lb, 14.40.Rt, 13.20.Gd 



I. INTRODUCTION 



The charmonium hadronic and radiative decays into light hadrons have provided an important way to probe the light 
hadron structures. In particular, with high statistics of J/yj and if) 1 events produced in e + e~ annihilation, the light 
Oh hadron spectra can be studied closely and dynamic information concerning the light hadron properties can be extracted 
from their production and decays. During the past few years, there have been several new resonance structures with 
J PC = h observed by BESII and BESIII in J/ip and ip' decays. They could be candidates of radial excitation 
states of the pscudoscalar mesons 77 and 7/', or exotic states such as glucball, multiquark state or hadronic molecule. 
, For instance, the BES-II Collaboration first reported a resonance structure in J/i/> -> 7A"(1835) — > jifn + TT~ [1], 
which was later confirmed by the BESIII measurement [2] with high statistics. Nevertheless, two additional resonance 
structures were identified as X(2120) and A"(2370) in the 77/7^ invariant mass spectrum [2, 3]. 

In fact, our understanding of the isoscalar spectrum is still far from well-established. Historically, the study of 
I the nature of r/(1405) and r/(1475) has been a hot topic and closely related to the effort of searching for the ground 
state pseudoscalar glueball in experiment. Since the first radial excitation states of 77 and 77' are generally assigned to 
7/(1295) and 7/(1475) taking into account their production and decay properties [4], it leaves out the abundant 7/(1405) 
. as a possible candidate for the pseudoscalar glueball. However, we would like to emphasize that such an arrangement 
still needs further studies, and it is still controversial whether 7/(1405) and 7/(1475) are two separated states or just 
one state of h in different decay modes [5]. 

With the availability of high-statistic J/ip and ip' events from the BESIII Collaboration, it allows us to tackle 
the question on the nature of 7/(1405) and 7/(1475). One important experimental progress is that the BESIII 
Collaboration [6] report the observation of anomalously large isospin violations of the 7/(1405/1475) — >■ 3ir in 



J/tp —> 77/(1405/1475) — > 77r°/o(980) — >• 7 + 3tt which, however, can hardly be understood by treating them as 
cither glucball or qq state. Interestingly, this decay process also explicitly involves the issue of ao(980)-/o(980) mix- 
ings. The BESIII data show that the Jo (980) signal is only about 10 MeV in width and the lineshape is different from 
the Breit-Wigner width of about 40 ~ 100 MeV [4]. Moreover, the isospin violation turns out to be significant with 
Bi?(r/(1405) -> /o(980)7T° -> 37r)/.Bi?(?7(1405) -> ag(980)7r° -> t/tttt) = (17.9 ± 4.2)%, which cannot be explained by 
the ao — /o mixing intensity measured in other channels [7] . 

An immediate theoretical interpretation is given by Rcf. [8] , where we propose that a triangle singularity mechanism 
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(TSM) via the intermediate K*K + c.c. rescatterings would lead to significant enhancement of the isospin violating 
decay, i.e. 77(1405/1475) —> K*K + c.c. —> /o(980)7r. In this transition, the dominant contributions would come 
from such a specific kinematic region that all the intermediate mesons in the triangle loop are literally on-shcll. The 
identification of such a mechanism seems to be nontrivial since it can naturally explain the narrow width of the /o (980) 
observed in the two pion invariant mass spectrum. Consequently, it raises an essential issue concerning the nature 
of 77(1405/1475) since the TSM can also contribute to the decays of 77(1405/1475) — > KKtt and r/Trir, and distort 
the lineshapes and shift the peak positions of the ?7(1405/1475) in those decay channels. As a result, a coherent 
study of 77(1405/1475) — > KKtt, r/irir, and 3ir is necessary and could be a key towards a better understanding of the 
7?(1405/1475) puzzle. 

In this work, we shall provide a detailed analysis of 77(1405/1475) — > KKtt, rjirir, and 3ir. We shall show that only 
one h isoscalar state, namely 77(1440), is needed in this mass region. With this "one state" assumption, we shall 
demonstrate that the TSM can lead to different mass spectra for 77(1440) — > KK* + c.c, ao(980)7r°, and /o(980)7r. In 
J/ip -> 777(1440), with 77(1440) -> KK* + c.c, a (980)7r°, and / O (980)7r, another possible contribution to the same 
final states is via /i(1420). Since the mass of /i(1420) is similar to that of 77(1440), we should investigate the role 
played by /i(1420) in these processes. In particular, due to the similar masses between /i(1420) and ?7(1440), the 
decay of /i(1420) would also experience the TSM. Therefore, a helicity analysis of the invariant mass spectrum for the 
overlapping /i(1420) and 7/(1440) is necessary. In comparison with the results reported in Rcf. [8], we have detailed all 
the analysis by including the /i(1420) contributions. Wc confirm the BESIII results by detailed helicity analysis from 
which we can extract the invariant mass spectra for 77(1440) in different channels. These features as a consequence of 
the TSM could be a natural solution for the long-standing puzzle about the nature of ?7(1405/1475) in experimental 
analyses. 

We also mention that the 77(1405/1475) — > 3ir decay was also studied in Ref. [9] recently in a chiral unitary approach. 
By exhausting several models and taking constraints from the meson-meson scatterings, the authors confirm that only 
ao(980)-/o(980) mixing can not explain the BES result [6] and the inclusion of the triangular diagrams is necessary [8]. 

The rest part of this paper is organized as follows: the formalism is presented in Sec. II. Section III is devoted to 
the numerical results and discussions. Our conclusion is given in Sec. IV. 



II. FORMALISM 



A. Effective Lagrangians and transition amplitudes 

The effective Lagrangians for the 77(1440) production have been presented in Ref. [8]. Here, we include the /i(1420) 
contribution and list the effective Lagrangians as the following: 

£viV 2 p = 9v 1 v 2 p£p.v P oPv 1 Pv 2 ^v 1 l l ) v 2 '>Pp , (!) 

£vp 1 p 2 = 9vp 1 p 2 (ippid fl ipp 2 - ippzdfti/jpjipfr , (2) 

£spiP 2 = 9sp 1 p 2 ipsi>p 1 ^p 2 , (3) 

£avp = gAVP^A^V^P , (4) 

= gie^ pa d»r^ p i r h +92e^ P ad^d x d»V i r h , (5) 

where S, P, V and A stand for four types of fields: scalar, pseudoscalar, vector and axialvector, respectively. For 
77(1440) the same diagrams as in Ref. [8] are calculated, while for /i(1420) the similar diagrams are listed in Fig. 1. 
Figure l(la)-(lb) are for /i(1420) KKn through K*K and a (980)7r channels. Figure l(2a)-(2b) are for /i(1420) -> 
7T + 7r~7r° through the TSM and ao — /o mixing. Figure 1(3) is for /i(1420) — > r]ir°ir°, where we assume that aQ(980)7r° 
gives the main contribution. 

Our kinematics conventions are shown in Fig. 2. Some common functions are defined as follows: 

Gf = s-m}+i^T f { S ) ' (6) 
G a = o -■ , ; , (') 
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FIG. 2: Kinematics defined in our formalism. 
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Then, the typical loop integrals can be expressed as 
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(11) 
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Taking into account that the relative signs between the charged and neutral loops are positive in isospin-conserving 
processes but negative in isospin- violating processes, it is convenient to define 



C p9k*Kir + C p9k*K-n > 

r c n c — r n n n 

c p9k*K-K c p9k*Kir > 

c pi5fc*if7r + c p!9k*K-K 

C p 1 9k*Kir ~ C p 1 9k*Kir 



(14) 
(15) 
(16) 
(17) 



where the superscripts 'c' and 'n' denote the charged and neutral loops, respectively. We also define 7^ , and dp 2 

for If it, in a similar way. 

The invariant amplitudes in Fig. 1 can then be expressed as 
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25/x/?5/„|G / ( S23 )| 2 |Zoop2 c - loopT 
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Mi = gie^ P y^ i+ g 2 e^py fi eP^ ■ Pfl B 2 (Q) . (28) 

In Eq. (27), we have put the identical factor 1/2 in the squared amplitude. The parametrization of Eq. (28) is taken 
from Ref. [10] and B 2 (Q) is the Blatt-Weisskopf barrier factor 

MQ) = ^Q* + WQl + 9Q* ' (29) 

where Q is the decay momentum, and Qo is a hadron scale parameter Qq = 0.197321/i? GeV with R the radius of 
the centrifugal barrier in fermi. In this paper we adopt R = 0.35 fm which is about the radius of J/ip. 



1. Helicity amplitudes 

In experiment, the quantum number of an intermediate state X is generally determined by measuring the angular 
distribution of the X decays. To proceed, we first make a model-independent analysis of the helicity structure of the 
transition matrix element which would allow us to separate different partial waves. Then, by comparing with the 
angular distributions measured by experiment, we can extract the dynamic coupling strengths for different partial 
waves. 

For a decay process a — > b + c with spin, helicity and parity (sj, Ai, r)i)i =a ,b,c, the decay amplitude in the rest frame 
of a can be expressed as [11], 

JW£ Ac (M;A a ) « D s x: * Xb _ Xc (ct>,e,0)F° :Xc , (30) 

where D x i * x x (0, ^, 0) is the rotation function, and F x a x is the helicity-coupling amplitude which is independent 
of angular variables. It satisfies two constraints taking into account angular momentum conservation and Parity 
conversation: 

|A 6 - A c | < s a , (31) 
nX = r) M (-y«- Sb - s °FX_ K . (32) 

From these relations, we can find out the independent helicity-coupling amplitudes. 

As follows, we will re-analyze the angular distributions of the recoiled photon in J/ip — >• -fX and the recoiled /o(980) 
in X —> /o(980)7r decays which are measured by the BESIII experiment [6]. Now as a preparation, we derive the 
helicity amplitudes of each vertex from the Lagrangians in Eqs. (l)-(5). 

In J/ip — > 7/1 , the helicity amplitude with transversely polarized /1 can be expressed as 

(A 7 ,A /l =±1|S'|A^) = gie^paPfc^e™ + g2B 2 (Q)e^ P aP^p"£* 1 p e}°ei, ■ Pj 

= -i 5l A 7 77ty£>i; o (<£,0,O) + ig 2 B 2 (Q)X 1 m 4 ,Q 2 D{l (cj ) ,e,0) 

= D{l (j>,9,0)Fl jXfi , (33) 

which is nonvanishing with A 7 = A f 1 , and 

F \-,x fl = -?'.9iA 7 m^ + ig 2 X^m^,Q 2 B 2 (Q) . (34) 

So the independent amplitude is 

Fft = im^-gx + g 2 Q 2 B 2 (Q)] . (35) 

When fi is longitudinally polarized, the g 2 term will have no contribution. The helicity amplitude is 

(A 7 , X fl = 0\S\X 4 ,} = gie^ upa p^e^e* p e* f l 

Et 

= -ig-LXym^, — -Dl* x U,e,0) 



with 



Dl: x (<j>,e, 0)Fl o , (36) 



Et 

F x o = -i.giA 7 m V; — - . (37) 
m fl 
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So the independent amplitude is 



Fw = -igim^^- . (38) 
For /i — > /o(980)7T°, the helicity amplitude can be written as 

(fon\S\fl) = 2ghK*Kgf KK 9 C K>K7rIfl ~ 9K'KttI}\ 

= 2g flK * K g foKK e%(p f _ l c~ +p lfl c ] ; i ) 

= D{* i0 (cf>,6,0)F$ , (39) 
where we have used the relation e p ■ p = 0. Also, in the helicity frame of /i(1420), we have 

e p -p 1 = -Dl* (<t>,9,0)Q . (40) 

So the independent amplitude is 

^00 = — %9fi K* Kg f a KKC^Q ■ (41) 

Similarly, the helicity amplitude for J/ip — > 77/(1440) can be obtained with both J/ip and 7 transversely polarized: 



(A 7 r?(1440)|5|Av,) = s^^fp^J?^^. 



fi n v c p ^er 
7 
1* 



= -iA^^m^X^*^ (0,0,0) 
= D\; Xy (cf > ,e,0)Ft rO , (42) 

with 

F\°o = -iX-yg^mm^ql . (43) 
For 77(1440) -> /o(980)tt , the helicity amplitude is 

(/ 7r|5|r ? (1440)) 

= 2gmK*K9f KK(gK* K7r Ifi ~ gK-K^l) 

= . (44) 
B. Angular distribution 

By combining the two-body decay amplitudes in the helicity frame, we can derive the total helicity amplitudes for 
the chain process J/ip — > jX — > y/o(980)7t° as shown in Fig. 3 and extract the angular distributions to compare with 
the experimental data [6]. 



. /o(980) 



FIG. 3: The kinematics for the chain process J/ip — > — > 7/o(980)7r°. 



For X being 77(1440), the total hclicity amplitude can be expressed as 
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The angular distribution is 
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where dfi = dr2 7 dil/ with df2 7 = dcos# 7 d0 7 and df2y = dcos#/ d</>/ . One step further, we obtain 

-^i- oc l + cos 2 7 , 
d cos o~ 
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dcos 6f 



oc const , 



(46) 

(47) 
(48) 



for the angular distributions of # 7 and #/ (98o)i respectively. These expressions are the same as those adopted in 
Ref. [6]. 

For X being /i(1420), the total helicity amplitude of the chain process can be expressed as 
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s - +im fi r fl 



lb 
00 



(49) 



The angular distribution via /i(1420) is 
dcr /i ~ 
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(50) 



where a = aie l ^ a is the ratio of the Xf t = ±1 amplitude to that of Xf 1 = 0. By integrating over corresponding polar 
angles in the above double distribution, one has access to the angular distributions of # 7 and 9 f , respectively, for the 
intermediate /i(1420): 
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(51) 
(52) 



With both 77(1440) and /i(1420) contributing to the chain process, the total helicity amplitude can be obtained in 
a similar way, namely, 
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where we have applied the selection rule A ft = A 7 for the transversely polarized /i(1420); r = r\e % ^ r is the ratio of 
the 77(1440) amplitude to that of /i(1420) with Xf 1 = 0. Then the angular distribution becomes 



da 



Vl+fl 

dO 



oc \ A m+h 

(A^,, A 7 )| 

= r 2 (cos 29 1 + 3) + 2r\ct\ sin#/ sin2# 7 cos^/ cos(0 Q — <f) r ) + 2r\ cos9f cos0 r (cos2# r + 3) 

+2a\ sin 2 9f sin 2 7 + ct\ cos0 Q sin20/ o sin2# 7 cos</>/ + cos 2 #/ (cos 29 ^ + 3) . (54) 

It is easy to show that Eqs. (46) and (50) can be reproduced by setting the corresponding resonance couplings to 
vanish. Similarly, the angular distributions of # 7 and 9f can be obtained 
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Checking Eqs. (47), (51) and (55), one can see that the cos6* 7 distribution is always symmetric as a feature of a two- 
body decay. In contrast, the angular distribution of cos#/ turns out to be nontrivial. As shown by Eqs. (48), (52) and 
(56), the contributions from different states with different quantum numbers are encoded in the angular distribution 
of cos9f . By fitting the experimental data, the coupling parameters can thus be determined which alternatively 
would provide information about the contributing resonances. As shown by Fig. 3 of Ref. [6], the cos#^ distribution 
is apparently asymmetric which indicates some contributions from the /i(1420) production besides ^(1440). 
With the explicit total helicity amplitude in Eq. (53), we can express the differential width as 



dF = (2^)5 16^J \P7\\Pfo\ d V^ d ^ d ^fo ■ 

We can define the following quantities by integrating the invariant mass yfsx of the /o(980)7r°: 
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which can thus be related to the quantities measured in experiment, i.e. 
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(58) 
(59) 
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From these relations, we can extract the information about the couplings from the angular distribution analysis. 

To compare with the experimental measurement of the unpolarizcd partial decay width in terms of the recoiled 
energy s by the photon in J/ip — > 7A — > ~/ABC, the following standard expression is adopted, 



dr,//^^ 7 jf^ 7 Asc _ 2sTj/^^ 7 x(s) x T x ^>abc(s) 



(62) 



where s is the four-momentum square of X = 7y(1440)//i(1420) in the reaction. A constant width r(/i(1420)) = 
0.0549 GeV is adopted for /i(1420) [4], while for 7/(1440), both constant width and energy-dependent form are adopted, 



1^,(1440) (s) 

where Ti a corresponds to Fig. l(la). 



r t (1440)^K*K^KK 
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III. RESULTS AND DISCUSSIONS 



In this part, we present our analyses and numerical results. First we demonstrate explicitly that the TSM is 
dominant in 7/(1440) — > 3n, and the main contribution is indeed from such a kinematic region that all the internal 
particles are close to their mass shells. Then, we give the fitting results about 7/(1440) from KKtt spectrum and 
show the predictions of ir + ir~ir° and r/ir ir a which are consistent with Rcf. [8]. By including /i(1420), we extract 
the couplings of /i(1420) through the analysis of the angular distribution of 7r + 7r~7r° channel [6]. Finally, we show 
that the combined results for both ry(1440) and /i(1420) in comparison with the BES data would allow us to draw 
a conclusion on the anomalously large isospin violations observed in 7/(1405) — > Sir and the nature of 7/(1405) and 
7/(1475). 



A. Loop integral 

Here we discuss in detail the calculation of I v i and Ifi in Eqs. (11)-(12). What we actually need in Ifi is the 
coefficients c p and c v \. We use two methods to calculate the loops: 

1. We directly calculate I v i and Ifi by LoopTools without any form factors. The UV divergences are regularized 
dimensionally by A = 2/(4 — D) — -fE + log47r, where A can be adjusted and the default value is A = 0. 

2. The exponential form factor method as applied in Rcf. [12]. Namely, an exponential form factor as follows is 
included in I v i and Ifi to cut off the UV divergence: 



exp 



k 2 — m 2 {p — k) 2 — m 2 (k — pi) 2 — m 2 



A 2 A 2 A 2 



(64) 



where A is the cutoff energy and characterize the effective range of the interaction. In principle, other forms of 
form factors can also be examined and we find the results are similar to each other. 

We present the calculations based on the above two treatments in Fig. 4. In the kinematic region that all the internal 
particles are close to their mass shells, namely the TS kinematics, these two treatments give nearly identical results 
since the form factor corrections are nearly unity. In particular, the absorptive part is dominated by the contributions 
from the TS kinematics. The real part turns out to be more sensitive to the form factor corrections when the internal 
particles deviate from their mass shells. Similar results are found for the /i(1420) since its mass is nearly the same 
as 7/(1440) and they share the same TSM. As a result, one can imagine that there should be little difference between 
these two treatments in the isospin- violating decay of t/(1440)//i(1420) — > 37r since the main contribution is from 
the absorptive part in the TS kinematics and the dispersive part would largely cancel out between the charged and 
neutral loop amplitudes. It is worth noting that the cancellation between the charged and neutral loop amplitudes 
eventually makes the calculation almost independent of the model uncertainties as explicitly pointed out in Rcf. [8] . 
This should be a direct way to confirm the dominance of the TSM in rj — > 3tt as a dynamic mechanism. 



B. Angular distribution analysis 

In the numerical calculations, the common coupling constants present in the triangle loops for 7/(1440) and /i(1420) 
are adopted the same as in Ref. [8], i.e. g a KR = 3-33 GeV, g am = 2.45 GeV, g K , K ^ = 3.208, and g K <. K ^ = 3.268. 
BES [13] and KLOE [14] give different values for the /o(980) coupling, namely, gj K x = 4.18 GeV and gf nn = 
1.66 GeV from BES [13]; and g fKR = 5.92 GeV and g /7r7r = 2.09 GeV from KLOE [14]. Similar to the treatment in 
Ref. [8], contributions from Figs. l(lb) and (2b) are neglected since they are only about 1/10 of Figs. l(la) and (2a), 
respectively. 

We adopt the mass and width of /i(1420) from the PDG [4], i.e. m A ( 142 o) = 1.4264 GeV and r /l(1420) = 54.9 MeV, 
but leave the mass and width of 7/(1440) to be fitted by the experimental data based on the "one-state" assumption. 
This is reasonable since the fi spectrum does not suffer from the ambiguity of possible abundant states in this energy 
region and as we shall see later that the /o(980) angular distribution measured by BESIII [6] only requires a small 
contribution from the /i(1420). 

Taking into account the present datum status, our analysis strategy is as follows: we first fit the BESIII data [6] 
for the angular distributions of the recoiled photon and / (980) in the decay of J/tp — > -fX and X — > / (980)7r°. 
This allows us to extract the relative strengths between the 7y(1440) and /i(1420) in the isospin violating decays. 
Then, by applying for the constraint from the J/tp -> 77/(1405/1475) -> jKKn from DM2, MARK III, and BES [4], 
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FIG. 4: The imaginary and real parts of the loop integrals in terms of the invariant mass of KK at y/s — 1.42 GeV. There are 
two sets of lines identified by those two spikes which correspond to the charged (lower mass) and neutral (higher mass) KK 
thresholds. For each set of lines, the thick and thin solid ones represent the results with A = 1.0 and A = 0.5 GeV, respectively, 
while the dashed lines denote the results of LoopTools calculation without form factor. 



it allows us to determine the absolute differential widths for both 77(1440) and /i(1420). We shall compare this 
with the exclusive fit by 7/(1440) as shown in Ref. [8]. In the end, we shall output the invariant mass spectra for 
7/(1440) / fi (1420) — > KKtt, rjwTr and 3ir, from which we would expect to observe different lincshapes and peak positions 
from the same state in different decay channels. 

The B 1 and 8f angular distributions of exclusive 7/(1440) and /i(1420) have been analyzed in Ref. [6] as parame- 
terized in Eqs (47)-(48) and (51)-(52). Now we consider the combined contributions from both 7/(1440) and /i(1420) 
and fit the BESIII data [6] using Eqs. (55)- (56) which can be expressed as 

6 7 (l + ccos 2 6» 7 ) , (65) 

b f (1 + ci cos 9 f + c 2 cos 2 6f ) , (66) 



dN 
d cos 6j 

dN 
d cos f 
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where 6 7 and 6/ are the overall normalization factors and 

= 1 - 2a\ + 3rf = 4ri cos^ 2 - g| 

C ~l + 2a?±3r 2 ' Cl _ a? + 2r? ' ° 2 ~ a\ + 2r\ ' 1 j 

We use the CERN program MINUIT to fit the data and the fitting results are demonstrated in Fig. 5. To compare 
with the results of Rcf . [6] , we show the x 2 values of different fits in Table I. From Fig. 5 we can see that the angular 
distributions are improved significantly when both 77(1440) and /i(1420) are included. The apparently asymmetric 
behavior of the cos#/ distribution in Fig. 5 can be well explained as the interference between 77(1440) and /i(1420). 



300 




C0Sf V»> cose , 



FIG. 5: Fitting results for the cos#/ (980) distribution (left panel) and cos# 7 distribution (right panel), respectively. The solid 
lines are the results considering both 77(1440) and /i(1420), while the dashed and dotted lines are the results with exclusive 
77(1440) or /i(1420), respectively. 



TABLE I: The fitting qualities of different fits. 



immediate states 


X 2 /d.o.f for cos 


9 7 x 2 /d-o.f for cos8j 


77(1440) 


40.2/15 


26.8/14 


/i(1420) 


59.0/15 


26.4/13 


77(1440) and /i(1420) 


38.3/14 


19.8/12 



The fitted parameters are as follows: 

b y = 118.5 ± 8.8, c = 0.538 ± 0.312 , (68) 

from the # 7 distribution and 

b fo = 145.7 ± 10.7, ci = 0.314 ±0.128, c 2 = 0.141 ±0.317 , (69) 
for 0f o . By solving Eq. (67), we obtain 

a\ = 1.197 ± 1.090, n = 1.50 ± 0.89, 4> r = ±(1.27 ± 0.20) , (70) 

which will allow us to extract the coupling constants for /i(1420), i.e. gf 1 (u2o)K*K, 9u an d 9i- Coupling gf 1 (u2o)K'K 
can be directly obtained from the width of /i(1420) [4] by assuming that KKtt channel is dominant; 3^ 7r; (i44o)/ ' 9i ls 
related to the fitted ri; and gijg\ is related to the fitted a\. 

The fitted parameters of ?y(1440) are 771^(1440) = 1-42 GeV , r i? ( 1 4 40 ) = 67 McV. The extracted couplings are listed 
in Table II. There are two solutions for the ratio of the £>-wave coupling to S-w&ve coupling gilg\i i-e. the value 
—0.179 indicates the S'-wave dominant, while the value 0.970 indicates the D-wave dominant. The present precision 
of the experimental data seems impossible to distinguish these two solutions. From our fit we find that the ratio of 
/i(1420) to 77(1440) in the KKn channel is about 17.3%. 
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TABLE II: Couplings extracted from the angular distribution analysis. 

9»7(i440)/f*x 3.638 

Sj/^(i440) (GeV- 1 ) (1.59 ± 0.32) x 1(T 3 

g fl K» K (GeV) 2.282 ±0.054 

31 (4.4 ± 2.7) x 10~ 4 

g 2 / 9l (GeV- 2 ) -0.179i°-J°S or 0.970t°;£| 



r( J/^-»7/i->7ifif tt) 
T( J/i/i->7T)(144())-»7Kj<'7r) 



(17.3 ± 23.4)% or (17.1 ± 23.1)% 



With these couplings, we can predict the corresponding spectra and ratios for J/rfj — > 7/i(1420) with /i(1420) — > 
KKir, TT + ir~ir° and t]tt°tt as shown in Fig. 6 and Table III. The results obtained by those two values of (72/ffi are 
almost identical as demonstrated in Fig. 6(a) for the KKir channel. So, in other channels we only show the results 
with 52/31 = —0.179. The main features of spectra and ratios are similar to those of 7/(1440). 

In the 7T + 7r~7r° channel, the results of KLOE are larger than those of BESIII by a factor of about 1.11 due to the 
difference of gfKii an d 9f-mr extracted from these two experiments as mentioned earlier. Meanwhile, it shows that 
the partial width (or branching ratio) is insensitive to the form factor cut-off energy. This feature has been discussed 
earlier and it is because that the model uncertainties will be largely constrained by the cancellation between the 
charged and neutral loop amplitudes. Nevertheless, the dominant contributions to the isospin-violating decays are 
from the TS kinematics where the form factor effects are rather small. 

In the 7/7r°7r° channel, the results are sensitive to the integration methods and cut-off energies due to the contributions 
from the dispersive part in the loop integrals. When varying the cut-off A from 1.0 GeV to 0.5 GeV, the results change 
about 31%. 





FIG. 6: Predictions for the spectra dY(J/ijj — > 7/i(1420) — > jABC)/dy r s. Figure (a) is for the KKir channel where the solid 
and dashed line denote results with different values of gijg\- Figures (b) and (c) are for the 7r + 7r~7r° channel with gf KK and 
gf Q1Tn determined by the BES and KLOE data, respectively, and (d) for the r/7r 7r° channel. The thick and thin solid lines in 
Figures (b-d) correspond to the results with A = 1.0 and 0.5 GeV, respectively, while the dashed lines denote the results by 
the LoopTools calculation without form factor. 

As we expect that the TSM appears to be more significant in /i(1420) than in 7/(1440) since the coupling /i(1420) — > 
K*K + c.c. is in a relative S wave, while 7/(1440) — > K*K + c.c. is in a P wave. Taking the results with A = 1.0 GeV 
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TABLE III: The extracted results for /i(1420), where Yabc = r(J/V> -> 7/1 (1420) -> 771BC) and 7?asc = Labc/iVs 



coupling gg/gi = -0.179 GeV~ 2 gg/gi = 0.970 GeV~ 2 



channel T(keV) 




R 


T(keV) 




R 


^i?yr 2.67 x 10" 


-2 


1 


2.63 x 10" 


-2 


1 


LoopTool 3.18 x 10" 


-4 


1.19% 


3.18 x 10" 


-4 


1.21% 


7r+7r-7r°(BES) A = 1.0 GeV 3.06 x 10" 


-4 


1.15% 


3.06 x 10" 


-4 


1.16% 


A = 0.5 GeV 3.11 x 10" 


-4 


1.17% 


3.11 x 10" 


-4 


1.18% 


LoopTool 3.54 x 10" 


-4 


1.33% 


3.54 x 10" 


-4 


1.34% 


7r+7r-7r°(KLOE) A = 1.0 GeV 3.40 x 10" 


-4 


1.27% 


3.89 x 10" 


-4 


1.29% 


A = 0.5 GeV 3.45 x 10" 


-4 


1.29% 


3.45 x 10" 


-4 


1.31% 


LoopTool 6.78 x 10" 


-3 


25.4% 


6.74 x 10" 


•3 


25.6% 


777T 7r° A = 1.0 GeV 4.68 x 10" 


-3 


17.5% 


4.68 x 10" 


-3 


17.8% 


A = 0.5 GeV 3.24 x 10" 


-3 


12.1% 


3.24 x 10" 


■3 


12.3% 



as an example, the ratio of 7r + 7r~7r° to KKtt is 1.27% in /i(1420), while in 7i(1440) the ratio is 0.762%. The ratio of 
?77r 7r to KKtt is 17.5% in /i(1420), while in ?7(1440) the ratio is 6.61%. The contributions from the /i(1420) also 
affects the peak position as demonstrated in the next Subsection. 

C. Invariant mass spectra including both 77(1440) and /i(1420) 

With the parameters fixed as the above, we compare the spectra and ratios with experiment [6] in Fig. 7 and 
Table IV where both 7/(1440) and /i(1420) are included. The main features are consistent with Ref. [8] where only 
?7(1440) was considered. 

It_shows that the contribution of /i(1420) is much smaller than that of 77(1440) in J/ip 7?7(1440)//i(1420) 
•yKKir. However, one should be reminded that this is largely due to the suppressed coupling for J/ip 7/1 (1420). 
In contrast, the contribution from the /i(1420) is relatively enhanced in the r]Tr°TT channel than in KKtt because 
of the TSM. The most interesting scenario is that the lineshapes of the invariant mass spectra for the KKtt, t]tt tt 
and Tr°ir + 7r~ decays are very different from each other due to the presence of the TSM in the last two processes. 
Also, the interferences of the TSM have led to the shifts of peak positions in those three channels which describe 
the experimental data consistently. Such a phenomenon retains even with contributions from ?7(1440) exclusively as 
found in Ref. [8]. 

For the isospin- violating channel of ?7(1440)//i(1420) 3ir, the observation of the narrow /o(980) in the 7T7r 
spectrum can be regarded as a signature of the TSM. As being shown in Ref. [8] , the narrow peak is located between 
the charged and neutral KK thresholds as a residual contribution due to the isospin violation. The mass difference 
between the charged and neutral kaons gives rise to the nonvanishing amplitudes between the K + K~ and K°K° 
thresholds which has been a crucial mechanism for the ao(980) and /o(980) mixing. Beyond this scenario, what we show 
here and in Ref. [8] is that the TSM can further dominantly enhance the /o(980) production in ?7(1440)/,/i(1420) — > Sir 
which eventually explains the anomalously large isospin violations. 

As listed in Table IV, we can see that the ratios of ir + ir~ir° and r]ir ir to KKtt agree well with experiment. 
Meanwhile, one also notices that the relative contributions from /i(1420) to 77(1440) are quite different in different 
channels as we discussed in Fig. 7. Namely, the relative strength of /i(1420) to 77(1440) turns out to be more significant 
in the r]ir ir channel than in the KKtt. It is because the S- wave coupling of K*K to /i(1420) would allow a relatively 
enhanced contributions from the TSM in J ftp 7/1 (1420) -> 7?77r 7r than 77(1440). 

In Table V we present the ?7(1440) and /i(1420) peak positions extracted in those three decay channels. Due to the 
contributions from the TSM, the peak positions are shifted differently. It shows that the exclusive results for 77(1440) 
and /i(1420) respectively or the results with their combined contributions have a similar feature. Namely, the largest 
peak mass can be seen in the KKtt channel, while the smallest one is the 37T channel. This qualitative pattern fits 
well the experimental observations in these three decay channels. 

Because of the TSM, the peak positions in both rjirir and 37r channels would move towards the K* K + c.c. threshold 
which is about 1.39 GeV. The more significant the TSM contribution is, the larger the peak position shift would 
be. As a result of the TSM dominance in the ir + ir~ir° channel, the peak position observed in the 7r + 7r~7r channel 
has a lower value than that in the t]tt 7t channel. The importance of the TSM suggests that a partial wave analysis 
including the TSM is necessary. Such a mechanism may also have significant interferences with the background. As 
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FIG. 7: The spectra dT(J/ip -> 7X -> yABC)/dy/s including both 7/(1440) and /i(1420). In the ifiS'Tr channel, we also show 
the experimental data, i.e. the solid triangles, solid circles, hollow circles, solid pentacles and hollow pentacles from MARK 
III(ifgif ± 7r =F ) [15], BES(JCg J ftr ± 7r =F ) [16], BES(K ± K^tt°) [17], DM2(K^K ± tv 11 ) [18], and DM2(K + K~ tt°) [18], respectively. 
The thin line in Fig. (a) is the background. In the -k + tt~ty° and n + TY~ channels, we show the results given by the LoopTools 
calculation with gf KK and gf 0W7r determined by BES. In the 7/7r 7r° channel, we choose the results with A = 1.0 GeV. 



TABLE IV: The combined results for Rabc = T(J/ip -47X-1 <yABC)/T(J/4i -> 7X -> jKKtt) including both 7/(1440) and 
/i(1420). The experimental data [6, 16, 19, 20] are also listed. For /i(1420) we adopt g 2 /gi = -0.179. 



R 


channel Theory 


Expt. 


KKit 1 


1 


LoopTool 0.781% 




7T + 7r-7r°(BES) A = 1.0 GeV 0.746% 


(0.90 ± 0.39)% 


A = 0.5 GeV 0.752% 




LoopTool 0.878% 




7r+7r-7r (, (KLOE) A = 1.0 GeV 0.837% 


(0.90 ± 0.39)% 


A = 0.5 GeV 0.843% 




LoopTool 23.7% 




t/ttV A = 1.0 GeV 8.22% 


(7.8 ±4.6)% 


A = 0.5 GeV 5.28% 





a consequence, it will lead to different lincshapes for the 7/(1440) in different production channels. For instance, the 
peak position of the 7/(1440) in J/ip — > jrjirir is slightly different from that in J/t/j — > uitjittt [21]. Note that the results 
of Ref. [21] are given by simple Breit-Wigner fit instead of partial wave analysis. Further detailed analysis of this 
channel using partial wave analysis should include the TSM as an important underlying dynamics in order to extract 
the correct pole position for the 7/(1440). 

With m^( 14 4o) = 1.42 GeV and mj 1 ( 1 4 2 g) = 1.4264 GeV, we find that the peak position shifts in the /i(1420) decays 
are larger than in 7/(1440). The reason again is because of the relative 5"-wave coupling for /i(1420) — > K*K + c.c. 
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When combining 77(1440) and /i(1420) together, the largest peak position shift that we can achieve is about 24 McV, 
which supports our one-state assumption. Namely, the 77(1405) and t/(1475) may be just one state in different channels. 



TABLE V: Peak positions in different channels. 



peak position (GeV) 


KK-ir 


7V + 7T TV 




77(1440) 


1.433 


1.416 


1.426 


/i(1420) 


1.431 


1.411 


1.422 


77(1440) +/i(1420) 


1.432 


1.415 


1.425 



D. Radiative decays of 77(1405/1475) 

Our proposal that 77(1405) and 77(1475) are the same state would have an explicit consequence in the description of 
the radiative decays of 77(1405/1475) — ^ jV, where V stands for the light vector mesons <f>, p° and uj. In the one-state 
assumption, ?7(1440) would be the SU(3) flavor partner of ?7(1295) as the first radial excitation states of 77 and 77'. In 
Ref. [5], it. was commented that by assigning the 77(1475) to the SU(3) partner of 77(1295) as the ss dominant state would 
not be able to explain why the observed branching ratios BR(n(lA75) — > 7/0°)) is larger than BR(rj(lA7b) — > j4>))- 
Also, it was commented that the observation that the much stronger production rate of J/tp — > 777(1405/1475) than 
J/tp — >■ 77/(1295) seemed not be obvious taking into account the above question. However, in this Subsection, we shall 
show that the experimental observations can be sclf-consistently understood by treating 7/(1440) and 7/(1295) as the 
SU(3) flavor partners. This can be explicitly demonstrated as the following: 

By assigning 7/(1295) and 7/(1440) as the first radial excitation of 7/ and 77', we can organize them as the following 
mixtures between 7777 = (uu + dd)/^/2 and ss: 

?7(1295) = cos a?777 — sin ass 

77(1440) = sin 0:7777 + cos ass , (71) 

where a is the mixing angle. 

In the J/tp radiative decays, it is a good approximation that the photon is radiated by the charm (anti-)quark, and 
the light qq of h is produced by the gluon radiation. By defining the production strength for the qq of h as the 
following: 

g = (qq\H\J/^,j) , (72) 
one can express the production amplitudes for 7/(1295) and 7/(1440) as 

M(?7(1295)) = (v^cosa- Rsma)g , 

M(rj(U40)) = (v^sina + iJcosa)^ , (73) 



where R = (ss\H\J/tp,"f)/go is an SU(3) flavor symmetry breaking factor, and one can simply set it to be unity as a 
leading approximation. It can be easily seen that a proper value for the mixing angle a in the first quadrant would 
lead to a much suppressed b.r. for J/tp —> 77/(1295) than for J/tp —> 77/(1440). The value of a can be determined 
by the b.r.s measured for these two channels. For instance, if one requires that B.R.(J/tp — > 777(1440)) / B .R.(J / 'tp — > 
777(1295)) ~ 10, namely, the production of 7/(1440) is about one order of magnitude larger than ?7(1295), one would 
have 

B.R.(J/j> ->■ 77/(1440)) _ / g„(i440) \ 3 / \/2sina: + i?coscA 2 
B.i?.(J/7/;-> 77/(1295)) ~ U„(i295)J I ^ cos a -i? sin a J 



10 , (74) 



where (7,7(1440) and 9,7(1295) are three momenta of the pseudoscalars in the J/ip rest frame, respectively, with R = 1, 
one has a ~ 38°. Such a mixing scenario will have explicit predictions for the radiative decays of ?7(1440) — > 07, p 0r j 
and W7. 

Since <p and w are nearly ideally mixed to each other and p° has isospin-1, we adopt the flavor wavefunctions, <p = ss 
and w = nn, and p° = (uu — dd)/y/2. The ?7(1440) radiative decays are via Ml transitions where the quark spin will 
be flipped by the magnetic interaction. A standard operator in the quark model can be written as 



2J e i^i ■ ~^-y\<i>SXA) , (75) 
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where /li = e/2rrii is the magnetic moment of the ith quark, and \4>sXa) and \4>aXs) are the flavor-spin wavefunctions 
for 77(1440) and vector meson, respectively. The subscriptions S and A means that the corresponding wavefunctions 
are symmetric or anti-symmetric under the exchange of the first and second quark (anti-quark). The flavor and spin 
wavefunctions are defined in a standard way as the following: 

<j>s{ss) = (ss + ss)/V2, 
<j)s(nn) = (nn + nn)/V2 , 

XA = (tWt)/V2 (76) 

for the pseudoscalar state, and 

<j>A{4>) = (ss- ss)/V2 , 
<j} A {p°) = ((uu-uu)-(dd-dd))/2 , 
4>a{w) = ((uu — uu) + (dd — dd))/2 , 

xs = n, u, (n + ;t)/V2 , (77) 

for the vectors. 

One can easily work out the flavor-spin couplings for those three channels as follows: 

, e 

h <t>l = ~n C0Sa > 

3m, 



h o~ = sm a 

2m q 



h u -y = sin a , (78) 

where m q = m u = and m s ~ 5m g /3. Apart from the spacial form factor and phase space factor in a P wave, the 
b.r. fraction among these decay channels are then 

2 -2 -2 

„ „ , „ „ , „„/ s cos q sm a sm a 

B.R.(-f4>) : B.R.( 1P °) : BJLfru) ~ : ^ r : . (79) 

For a proper value of a in the first quadrant, the decay of ?7(1440) — > jp° would be dominant. To be consistent with 
the production of 77(1440) and 77(1295) in the J/ip radiative decays, i.e. a ~ 38°, one obtains B.R.(j4>) : B.R.(jp°) : 
B.R.^uj) ~ 1 : 3.8 : 0.42. 

In brief, given a proper mixing angle between the 77(1295) and 77(1440) as the first radial excitation states of 77 and 
77', the theoretical interpretation of the 77(1405/1475) as a single state of 77(1440) does not obviously conflict with the 
so far available experimental data at all. The misunderstanding that the branching ratio of 77(1475) — > 70 should 
be larger than that of 77(1475) — > 7/5 if ?y(1475) is the higher mass partner of 77(1295) is not necessary at all due to 
the suppression of the quark masses in the Ml transition. This point, unfortunately, has not been realized in earlier 
analyses. 



IV. SUMMARY 



In summary, we have made a systematic analysis of the correlated processes J/ip —> 777(1440)//]. (1420) with 
?7(1440)//i (1420) — > KKw, rjTTW and 37r, where the role played by the TSM is clarified. Our combined analysis 
including 77(1440) and /i(1420) agrees well with the experiment data, and provides an overall description of the 
processes J/ip — > jX with X — > KKir, r]ir°ir , and 7r + 7r~7r°. In particular, we show that the inclusion of the 
/i(1420) can improve the description of the /o(980)7r° angular distribution significantly, although the contribution 
from /i(1420) is much smaller than 77(1440). By fitting the BESIII data for J/ip -> jX -> 7/ (980)7r°, we extract the 
coupling parameters of /i(1420). It allows us to estimate that the ratio of /i(1420) to ?7(1440) in the KKtt channel 
is about 17.3%. This does not change the results of the previous work [8] in which we assumed that 77(1440) was the 
only contributing state as treated by the BESIII. We also show that /i(1420) can contribute some percentages to the 
narrow peak of / (980) -» irir via the TSM. 

We emphasize that the dynamic feature of the TSM can be recognized by the strong narrow peak observed in the 
37T channel with the anomalously large isospin violations. Moreover, it leads to the obvious peak position shifts for 
the same ?7(1440) or /i(1420) state in different decay channels, which may suggest that the 77(1405) and 77(1475) are 
actually the same state. So far, such a one-state prescription seems not to have a conflict with existing experimental 
data. This may shed a light on the long-standing puzzling question on the nature of 77(1405) and 77(1475) in the 
literature. 
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